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Abstract 

We explore extensions of inelastic Dark Matter and Magnetic inelastic Dark Mat- 
ter where the WIMP can scatter to a tower of heavier states. We assume a WIMP 
mass m^^ ^ 0(1 — 100) GeV and a constant splitting between successive states 6 ~ 
0(1 — 100) keV. For the spin-independent scattering scenario we find that the direct 
experiments CDMS and XENON strongly constrain most of the DAMA/LIBRA pre- 
ferred parameter space, while for WIMPs that interact with nuclei via their magnetic 
moment a region of parameter space corresponding to ~ 11 GeV and 5 < lb keV 
is allowed by all the present direct detection constraints. 



1 Introduction 



The nature of dark matter is one of the fundamental questions facing physics, with numerous 
experiments being performed to search for it both directly and indirectly. Direct detection 
experiments hope to detect Weakly Interacting Massive Particles (WIMPs) by observing 
the recoil of target nuclei which interact with ambient WIMPs in the nearby galactic halo. 
One particular approach boils down to a counting experiment, looking for a signal in excess 
of all known background processes. Due to the expected small scattering cross-sections 
involved, these experiments typically need to reject many large backgrounds such as natural 
radioactivity and cosmic rays in order to be sensitive to dark matter scattering. 

Another approach to direct detection of dark matter relies on the relative motion of 
the Earth through the dark matter halo. As the Earth orbits the Sun, which in turn is 
moving about the center of the Milky Way, the flux of WIMPs impinging on the Earth 
undergoes an annual modulation [1]. In particular, it is more likely to find WIMPs at high 
relative speeds when the Earth's motion is aligned with the Sun's (in the summer) than 
when they are pointing in opposite directions. A larger relative speed results in an increased 
potential for higher energy nuclear recoils. Given the finite energy thresholds of direct 
detection experiments, the result translates directly into an annual modulation of the rate 
of WIMP scattering with target nuclei. This modulation helps isolate a WIMP signal from 
the known backgrounds, which are not expected to display a strong modulation, and allows 
for detection of WIMPs without having to identify individual events as arising from signal 
or background. In fact, the DAMA/LIBRA experiment has reported evidence for just such 
an annual modulation signal whose peak is in June, consistent with the expectations of dark 
matter scattering IJJ . However, in the simplest dark matter models the typical cross-sections 
needed to generate a modulation signal large enough so as to explain the DAMA/LIBRA 
results are so large that they are inconsistent with the null results at other direct detection 
experiments, including XENONIOO [3] and CDMS II [3]. 

Inelastic Dark Matter (iDM) [5i 16j models attempt to resolve this puzzle. iDM models 
alleviate direct detection experimental constraints by requiring inelastic scattering of the 
form 

X + N^X* + N (1) 

where x is the incoming WIMP, x* is a heavier state into which it must scatter when interact- 
ing with a target nucleus A^. Only scattering which transfers at least energy 6 = m-^^ — is 
kinematically allowed, which gives more relative weight to large velocity scattering, enhanc- 
ing the annual modulation effect. In addition, since the momentum transfer of the scattering 
is controlled in part by the target mass (assuming it is not much greater than m^), the scat- 
tering rate is also different for different mass target nuclei. Initially, for ~ 100 GeV and 
6 ~ 100 keV, iDM models were able to explain the DAMA/Libra observations, while not 
being ruled out by other experiments. In the time since they were proposed, the increases 
in sensitivity have since closed the window of parameter space which could explain DAM A. 
For example, XENONIOO has ruled out the iDM parameter space with a 5 ^ 120 keV that 
explains DAMA as scattering off of Iodine at the 90% confidence level [3 [HI [9] . 
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In addition to its atomic mass number, another feature which distinguishes target nuclei 
are their electromagnetic properties, including charge and magnetic moment. Magnetic in- 
elastic Dark Matter (MiDM) [T^ exploits these differences by positing an inelastic WIMP 
whose primary interaction portal with Standard Model (SM) particles is by exchanging pho- 
tons by virtue of a magnetic dipole moment [HI [121 1131 [lH [Ej . Such a photonic portal favors 
target nuclei with larger charge and/or magnetic moment, and leads to an enhanced rate at 
DAMA, as compared to CDMS (for example) since Sodium and Iodine have high magnetic 
dipole moments in comparison to Germanium (see Table [T|. A viable MiDM parameter space 
results y^|16j, subject to mild (and somewhat model-dependent) constraints from the null 
results of the Fermi/GLAST search for gamma-ray lines [IT] from WIMP annihilation 
and from LEP searches [T9| for missing momentum [121 EDI 121] • 

Aesthetically, the need for a small splitting S relative to other mass scales in the theory 
such as is somewhat mysterious. The existence of such a splitting may be motivated 
by introducing a weakly broken symmetry [TUl 12^ which would otherwise require the 
elements of a WIMP multiplet to be degenerate, or can be produced in models where the 
WIMP is a composite state, bound either by a confined [211 [25] or by a weak long range 
force p6|. Such models naturally accommodate multiple mass scales, but it remains true 
that one tunes a parameter in order to generate a fine or hyper-fine splitting of the correct 
size to generate a realistic iDM model. 

In this article, we explore a variation of composite inelastic models with a new feature 
ameliorating the need to tune any parameter related to 6. We consider a class of models 
in which the WIMP is a bound state of a new confined gauge force, where the dark matter 
is the lowest lying state consisting of a heavy preon which provides the bulk of the WIMP 
mass (and perhaps "flavor" quantum numbers which insure its stability), as well as some 
other light fundamental degrees of freedom. The confinement scale A of the new gauge sector 
satisfies A ^ m^, such that one can expect a continuum of excited states whose levels are 
much smaller than the WIMP mass itself. A is also much smaller than the typical scale 
required of iDM in order to explain the DAMA results. The collective up-scattering from 
ground state into a variety of the excited states will explain the DAMA modulation results, 
with the characteristic splitting scale emerging organically from a theory whose characteristic 
splittings are somewhat smaller. Since this dark matter candidate looks something like a 
massive preon surrounded by a rather flimsy (easily perturbed) gluosphere, we refer to it as 
"fluffy" dark matter (fDM). 

We assess the conditions under which fluffy dark matter can fit the DAMA modulation 
signal by interacting with nucleons either through a (canonical) Higgs-like portal, or by 
virtue of a magnetic moment through the photon portal. We model the fluffy spectrum as 
consisting of a ground state with a continuum of excited states at roughly evenly spaced 
intervals above it. As we will see shortly, the upshot is that the Higgs-like portal is ruled 
out by a combination of XENONIOO [7] and low threshold CDMS [27] data, but magnetic 
fluffy dark matter (MfDM) is viable for WIMP masses around ~ 11 GeV and 6 ~ keV. 
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2 Effective Theories of Fluffy Dark Matter 



Fluffy darlc matter can arise as the low energy limit of a variety of UV theories. While 
it would be interesting to pursue some detailed examples (allowing one to discuss colhder 
signals, which could show features of unparticles [28] or a hidden valley models [29]), we 
defer construction of detailed models to future work. Instead, we focus on the low energy 
dynamics, which we express in terms of an effective field theory containing the composite 
bound states. 

As an underlying picture, we imagine a large N confining gauge theory, with confinement 
scale A. We refer to the vector bosons associated with the new gauge sector as "gluons" 
(without confusion with respect to the force carriers of the SU{3)c of the SM). There are 
also matter fields consisting of a heavy adjoint Majorana fermion (the "gluino") which is 
a singlet under the SM gauge interactions, as well as some connector fields with both SM 
gauge quantum numbers, and transforming under the new force such that they can mediate 
interactions between bound states of the new gauge force and either the SM Higgs or photon 
by inducing a magnetic moment. It would be interesting to explore candidate models in 
terms of RS-like dual theories [30] , but we leave such constructions for future work. 

The lowest lying states of the confined sector consist of glueballs with masses of order A 
(which typically can decay into SM particles) and glueballinos, Xi, Majorana fermion bound 
states of gluinos and gluons. We assume an underlying flavor symmetry which renders the 
gluino, and thus the Xi^ stable. The WIMP is identified as the lightest of these states with 
mass m^, and there is a sector of excited states with masses larger by order A. We will assume 
that the spectrum of excited states consists of a tower of states labelled by an integers n 
with masses. 



where 5 ~ A ^ m^. 

The connector fields are responsible for inducing interactions between the bound state 
WIMPs and the SM fields. One such interaction is a scalar coupling to two Higgs doublets. 



where H is the SM Higgs. After electroweak breaking, this leads to an interaction of XiXj 
with a single Higgs boson with strength r^v. We have simplified the discussion by assuming 
that a single parameter controls the interaction, independent of i and j. In principle, this 
interaction can mediate elastic as well as inelastic scattering, but since we will find this case 
has difficulty explaining the DAMA results in the light of other null searches anyway, we will 
not dwell on this point here and just assume that something forbids the i = j terms from 
occurring. We refer to a model coupled in this way as fluffy dark matter (fDM). 
The second interaction we consider is a magnetic moment interaction. 



rrin = rriy + n 5 



(2) 



XiXj H^H, 



(3) 



(4) 
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where fi^ parameterizes the strength of the magnetic dipole, a^^ = 2[7^,7i,]/2, and F'^'^ is 
the photon field strength. For a Majorana fermion, this interaction vanishes for i = j, which 
nicely explains why the WIMP should scatter inelastically off of nuclei. We have again made 
the simplifying assumption that the same parameter /i^ at least approximately describes the 
strength of the interaction for all i and j. We refer to a fiuffy WIMP which interacts with 
the SM primarily in this way as magnetic fluffy dark matter (MfDM). 



3 Direct Detection of Fluffy Dark Matter 

3.1 Dark Fluffy Scattering 
3.1.1 fDM 

For an inelastic scattering, the differential scattering rate may be written [31 
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and Nt is the number of target nuclei per kilogram, p^^ is the local WIMP energy density, 
(j„ is the cross-section for scattering off a single nucleon, m^v is the nucleon mass, /i„ is the 
reduced mass of the nucleon- WIMP system, Er is the recoil energy, cos 7 is angle between 
the velocity of the earth and the recoil velocity of the nucleon as seen the earth's rest frame, 
F'^{Er) is the helm form factor describing the loss of coherence of the nucleus at large 
momentum transfer [32], vo — 220 km/s is the WIMP velocity dispersion, fesc — 500 km/s 
is the escape velocity, ve is the velocity of the earth, and n{vo,Vesc) normalizes the velocity 
distribution. We assume here a "standard" Maxwellian distribution of WIMP velocities in 
the halo. Conservation of energy and momentum dictate that the minimum velocity to 
scatter is, 
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In the case of fDM, the Wimp scattering is into one of the whole tower of excited states. 
For a given flnal state excited WIMP j one has, 
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where as discussed above, we take 5^ ~ jS and ~ cr„. This resuhs in three relevant pa- 
rameters controlhng the predictions: the WIMP mass m^, the sphtting between consecutive 
excited states 5, and the cross-section with nucleons (T„. 

3.1.2 MfDM 

In MiDM scenarios the interactions between the WIMP and the target nucleus are mediated 
by photons which couple to the WIMP's magnetic dipole moment The interaction with 
the target nucleus can either proceed through its charge or magnetic moment, leading to 
both dipole - dipole (DD) and dipole-charge (DZ) interactions. The differential cross-section 
with respect to the recoil energy is. 
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where is the heaviest state kinematically accessible [10]. A list of values for Z and Hn for 
a variety of relevant target nuclei is presented in Table [TJ The DD term is proportional to 
yU^^^, implying that the rate at DAMA (whose target is Nal) can be significantly enhanced 
as compared to experiments using target nuclei with lower /i, such as xenon or germanium. 

In Eqs. (10)-(12), F'^[Er] is the Helm form factor [32] and Fd^[Er] is the nuclear magnetic 
dipole form-factor. While the helm form factor is expected to be a good approximation for 
the factor by which the cross section decreases for non-zero momentum transfer in spin 
independent interactions, the nuclear magnetic dipole moment form-factor deserves more 
careful treatment. The usual thin-shell model for spin-dependent interactions need not be 
a good approximation for heavier nuclei [33]. We follow Ref. [3l] in obtaining the form 
factors for ^^^Xe, ^^^Xe, ^^^Cs and ^^^I. The theoretical model therein from which we have 
derived our form factor reproduces observables like magnetic moments of relevant nuclear 
states quite well. Nonetheless, as discussed in more detail in Ref. [TU], this remains a source 
of uncertainty in our predictions which is difficult to quantify in the absence of more direct 
experimental inputs. We define the form factor as 

F\Er) = ^ [flnlF^,{Er) + 2Uhn,n,F,,{E,) + fln\FUE,.)) (13) 
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Table 1: Natural isotopes with quantities relevant to direct detection searches [35] . 

where Qq, Qi are effective g factors (Table IV, Ref. [31]), and Fqq, Fqi, Fu are spin structure 
functions (Figs. 3-4, Ref. 01]). The coefficients fo = fin + l^p and /i = /ip — /in are the 
isoscalar and isovector couphng constants, jin = —l.dfiN and fip = 2.8/i7v are magnetic 
moments of the proton and neutron. The overall factor N is fixed by imposing the condition 
that the form factor is 1 at zero momentum transfer. 

For all other nuclei we use the distribution [35] which takes into account the coupling to 
all 'odd-group' nucleons , 

p2/ ^ N _ / Jiiir^) {qrn < 2.55, qr^ > 4.5) , . 

W^nj - I Q_Q4-^ ^2.55 < qrn < 4.5) " ^ ' 

3.2 Fits to the DAMA/LIBRA Signal 

Our next task is to assess the ability of fDM and MfDM to fit the DAMA observations. We 
scan over and 6, and construct a function to fit the amplitude over the 2 — 8 keVee 
energy bins reported by DAMA/LIBRA, using the quenching factors 0.3 (0.09) for Na (I) [37] . 
To construct the x^-function, we estimate the modulated amplitude as half of the difference 
between the maximal and minimal scattering rates during the yeai[^ We focus on the 
2 — 8 keVee region, since for higher energies both the observed and predicted spectra are 
consistent with zero modulation and the inclusion of these energy bins would act to artificially 
improve the quality of the fit. 

For each point in the m^ — S plane, we find the corresponding spin- independent scattering 
cross-section cr„ for fDM or the dipole moment fi^ for MfDM that leads to the minimal value 
of x^- The results are shown in Fig. [l| Using these values of (j„ for fDM and /i^ for MfDM 

^This approximation breaks down if the scattering rate were to be zero for a substantial part of the winter 
months. However we find that this approximation is vahd for aU of the DAMA/LIBRA preferred regions of 
parameter space. 
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Figure 1: Best fit values of cr° in units of 10~^° cm^ (left panel) in the fDM scenario and fi^ 
in units of lO'^fi^ (right panel) in the MfDM scenario. The maximum value in both cases 
is along the edge of the tan (off-white) region and is 24 in the fDM case and 2.8 in the fDM 
case. The white region in both plots is not allowed as the 6 is too high for any scattering to 
take place. The kink in the boundary of the allowed region at m-^ ~ 40 GeV is due to the 
fact that below this mass the upper limit on 6 is set by scattering off of Sodium, and above 
it the limit is set by scattering off of Iodine. 
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Figure 2: Constraints on fDM, for large {m^,5) (left) and small {m^,5) (right), from direct 
detection experiments. The darker blue (grey), blue (grey) and lighter blue (grey) regions 
correspond to the fDM points that fit the DAMA/LIBRA annual modulation spectrum at 
68%, 95%, and 99.7% C.L. In the left panel regions to the right and below of each solid 
black line are ruled out by the corresponding direct detection experiment. In the panel on 
the right, parameter space to the left of the line corresponding to the CDMS low threshold 
limit and to the right of the line from XENONIO '10 are excluded. 



we can calculate the scattering rates at a particular direct detection experiment using Eq.(|5]) 
and Eq.(lO) respectively. In this sense, the bounds on regions of and 6 that we later 
derive from other experiments are themselves dependent on DAMA. We consider bounds from 
XENONIOO (both the 11.17 days of data collected in 2009 [38] and the 100 days of data 
collected in 2010 [3]); CDMS II (including the low energy threshold analysis) [1]; KIMS [32] ; 
and COUPP [ID]. Bounds from XENONIO [41], ZEPLIN III [32], CRESST II [321 El], and 
PICASSO [35] were considered, but do not appear in the figures as they are weaker than the 
bounds from other experiments. A summary of the data used to constrain fDM and MfDM 
scenarios is presented in Table |2} 

In the plane of and S, we find islands where both fDM and MfDM can provide an 
acceptable fit within error bars to the DAMA/LIBRA data. We map out the contours of 68%, 
95%, and 99.7% C.L. agreement with the DAMA/LIBRA annual modulation amplitude, 
which we plot as the darkest blue (grey), dark blue (grey), and lighter blue (grey) shaded 
regions of Fig. [2] and Fig. [3] for both fDM and MfDM, and for low mass and high mass 
WIMPs. These contours are somewhat conservative in the sense that the best fit points 
have < 1 per degree of freedom. The tan (off-white) regions of each plot have poorer 
agreement, and the white regions of each figure have 6 large enough that fmin > Vesc, so no 
scattering is possible. 
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Figure 3: Constraints on MfDM, for large {m^, 6) (left) and small {tti^, 6) (right), from direct 
detection experiments. The darker blue (grey), blue (grey) and lighter blue (grey) regions 
correspond to the MfDM points that fit the DAMA/LIBRA annual modulation spectrum at 
68%, 95%, and 99.7% C.L. The parameter space below and to the right of each of the black 
lines are ruled out by the corresponding experiment. 

3.2.1 fDM 

In the case of fDM, the high WIMP mass region is shown in the left panel of Fig. |2} The 
regions that fit the DAMA/LIBRA data well are approximately the same as the those pre- 
viously identified in Ref. ^61 HTj as being preferred by iDM. As such, this parameter space 
represents a WIMP which is not particularly fiuffy because only one excited state is relevant 
for the scattering. It thus suffers the same fate as standard iDM and is highly constrained 
by XENONIOO and CDMS for a Maxwellian halo 

In the low mass region, there is a separate island of good fit region, where the scattering 
is dominantly off of Na nuclei. Here, we find roughly any 5 < 20 keV works to describe 
DAMA, which for small splitting looks very much like a fiuffy WIMP. However, WIMP 
masses above around 11 GeV are excluded by XENON 100 (which is essentially enough to 
remove the interesting parameter space) and the CDMS low threshold analysis removes the 
entire interesting low mass region as well. It is interesting to note that unlike the other 
exclusion limits, the low threshold CDMS limit excludes the region to its left, for smaller 
m^. This is due to the fact that the values of (T„ extracted from fitting to DAMA/LIBRA 
increase as decreases, as can be seen in the left panel of Fig. [T| which compensates the 
smaller typical scattering energy all the way to the region where no scattering is possible. 

Combining the low and high mass regions, we see that the regions of fDM parameter space 
consistent with the DAMA energy spectrum are strongly constrained by the null search limits 
from other direct detection experiments. 



9 



r ( Vn PTl m PTl f 
±UJ\.IJ CI llliCliu 


r jVnot;i 1 TP 


T^iTTiP PpTion 

J. llliC' -L C-llWU. 


kJlclliCLl VV llll.J.'^VV 


( ) nt;PTA/'pn 


UJJMo Uo 


iy4.1 Kgcl 


v/i /r\T n/i /r\Q 
i 1 i/U / - y/ 1/Uo 


iU - iUU KeV 


o 

L 


LjUNld Low-ii/nergy 


z4i KgQ 


1 m /C\P. O /I /HQ 

lUi/Uo - y/ i/Uo 


z - KeV 




XhjN (Ji\l lU 


U.o X 316.4 kgd 


10/6/06 - 2/14/07 


4.5 - 75 keV 


13 


AhjiN UiN iUU 


161 kgd 


1 n /on 1 c\c\ n /i o lc\c\ 

10/20/09 - 9/12/09 


T /I on 1 l..^"\7' 

(A - 29.1 keV 









^ fi /s /t n 

1/ LO j lU - \J 1 1 lU 


^ zL AA f> VpV 

O.^ - ^^.U JVC V 


Q 

o 


ZEPLIN III 


0.5 X 63.3 kgd 


2/27/08 - 5/20/08 


17.5 - 78.8 keV 


5 


CRESST II (W) 


0.59 X 0.9 X 48 kgd 


3/27/07 - 7/23/07 


12 - 100 keV 


7 


CRESST(O) 


564 kgd 


7/09 - 10/10 


~ 10 - 40 keV 


32 


CRESST(Ca,W,0) 


730 kgd 


7/09 - 3/11 


~ 10 - 40 keV 


67 


KIMS 


3409 kgd 


1/05 - 12/06 


20 - 100 keV 


955 


COUPP 


28.1 kgd 


11/19/09 - 12/18/09 


21- 200 keV 


3 


PICASSO 


13.75 kgd 


6/07 - 7/08 


20 - 200 keV 


2051 



Table 2: Exposure and data collection period, signal energy window, and number of observed 
events at various dark matter direct detection experiments. For CRESST we write ~ 10 keV 
for the sake of brevity in the table. The threshold is different for each of their detector 
modules. The Observed events at KIMS has been inferred using the exposure and event 
rates from Ref. 1391. 



3.2.2 MfDM 

The results for MfDM in the high WIMP mass range (m^ ~ 100 GeV) are shown in the left 
panel in Fig |3} The DAMA/LIBRA spectrum is well fit in the region m-^ ~ 80 GeV and 
b ~ 100 keV. The strongest constraints obtained are from XENONIOO and KIMS (whose 
Cs and I nuclei have relatively high nuclear magnetic moments, and are thus particularly 
good probes of MfDM). XENONIOO sets the stronger constraint and excludes the entire 
region allowed by DAMA. Similarly to the fDM scenario, this region of parameter space 
does not correspond to a very fluffy WIMP and the region of parameter space favored by 
DAMA/LIBRA is similar to that found in the MiDM scenario of Ref. [lO]. 

In the low WIMP mass region, we find that masses around 14 GeV and with 5 < 20 keV 
provide a good fit to DAMA. This region is shown in the right panel of Figure [3} where we 
see relevant constraints from XENONIOO and COUPP. CDMS low threshold constraint is 
weak, due to the small magnetic coupling, while the XENONIOO limit only constrains a part 
of the parameter space that provides a good fit to the DAMA/LIBRA data. The XENONIOO 
limit allows a significant portion of the MfDM parameter space that fits the DAMA/LIBRA 
signal at the > 95% C.L. The fit in this allowed region is also not dependent on 5 for values 
less than 15 keV. 

The larger 8 region indicates that MiDM with WIMP masses around 11 GeV is consistent 
with current direct detection data. The small 5 region realizes the hope of a fluffy WIMP, 
with many excited states contributing to the scattering. For example, an 11 GeV WIMP 
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can scatter up to states that are heavier by < 15 keV, so for a sphtting of (5 = 1 keV, 15 
excited states participate in the scattering. For smaller values of S, the number of relevant 
states in the scattering increases. 

4 Outlook 

We have considered a picture where dark matter is the lowest lying state in a composite 
sector, whose low energy scattering with nuclei is inelastic. We refer to such a WIMP as 
fluffy, since its internal structure as a confined state is rather easily perturbed by its environ- 
ment. For very low confinement scales, the preferred inelastic splitting emerges somewhat 
organically, driven largely by the energy thresholds of DAMA itself. We perform fits and find 
that light (~ 11 GeV) WIMPs with small splittings 5 provide the best fits to the DAMA ob- 
servation of annual modulation and energy spectrum. For a WIMP whose interactions with 
the SM are through iso-spin conserving spin-independent couplings to the SM, constraints 
(largely from XENON 100) are enough to close off the region of parameter space able to ex- 
plain the DAMA signal. However, a fluffy WIMP whose interactions are magnetic in nature 
can explain DAMA and still remain consistent with the bounds from other direct detection 
experiments. Because of the low masses favored by the fit, it is difficult for XENON to 
access the parameter space of interest. Lower threshold experiments perhaps offer the most 
promising probes in the future. 

A low compositeness scale is a challenge for cosmology, and leads to potentially many 
unusual features for a theory of dark matter. To begin with, our analysis assumed the bulk 
of the WIMPs in the halo were in the ground state. If this were not the case, the WIMPs 
would down-scatter as well as up-scattering. One could imagine engineering this possibility 
directly into a model of luminous dark matter [38], which is an interesting independent line 
of investigation beyond the scope of this work, but we prefer to imagine there are interactions 
which efficiently de-excite the WIMPs, perhaps through a coupling to neutrinos [49]. 

In the early Universe, a fluffy WIMP can undergo a freeze-out process which differs 
substantially from a standard WIMP. The seed partons may freeze out at early times, but 
the confined states may come back into equilibrium because of the surrounding clouds of light 
partons increase their effective cross sections after the phase transition. This can lead to an 
additional dilution of the WIMPs at late times [^01 ED E2]- Alternately, one could explore 
the case where excitations with large splittings continue to play an active cosmological role, 
perhaps obviating the need for an external stabilization symmetry [531 [53] • 

Since we favor a low confinement scale, there will be additional (at least) gauge degrees 
of freedom contributing to the thermal bath even at relatively late times. Precision mea- 
surements of the primordial abundances of light elements produced through big bang nucle- 
osynthesis, and of the cosmic microwave background, provide tight constraints on additional 
light degrees of freedom during the relevant epochs [551 EH] . Both sets of measurements are 
currently consistent with no new degrees of freedom being present, with error bars allowing 
for around at most one new state equivalent to an "effective neutrino species". Future CMB 
measurements by PLANCK are expected to have the precision to shrink the uncertainties 
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to the order of a tenth of an effective neutrino species One could imagine evading 

these bounds if the dark sector has a separate temperature from the SM plasma [52], or 
if the coupling of the new confined force is time-dependent, perhaps being set by value of 
a modulus through a term such as (pF'^'^F^^, where (0) starts at large values, leading to a 
tightly bound WIMP during the cosmologically sensitive times, but rolls to lower values in 
late cosmological times, leading to the WIMPs puffing up. 

After confinement, one can expect the analogues of glueball states for the new gauge 
force (whose masses should be roughly A ~ keV) could contribute to relevant phenomena. 
For example, if sufficiently long-lived, they could end up as a subdominant component of 
warm dark matter, or their decays could contribute additional entropy to the Universe. In 
addition, such states may have large (strong force residual) interactions with the WIMPs. If 
WIMPs can efficiently lose kinetic energy in collisions, either by converting kinetic energy to 
excitation energy and then de-exciting by radiating a glueball, or by exchanging glueballs in 
elastic collisions, it can cause elliptical galaxies to become spherical |SB1 [22] . One can even 
imagine more radical shifts in galactic dynamics, such as cases where some fraction of the 
binding is due to glueball exchange (in addition to gravity), with the galaxy itself looking 
something like the analogue of a heavy nucleus state of the new gauge force. 

One can imagine other scenarios in which a fluffy WIMP might provide an interesting 
model of dark matter. For example, models of exciting dark matter [60j invoke an MeV 
split excited WIMP to explain the INTEGRAL/SPI 511 keV gamma ray excess [61]. One 
could easily imagine a fluffy WIMP model allowing this small scale to emerge organically 
as it did here in an inelastic scattering context. It would be interesting to see if a common 
framework could potentially explain the DAMA signal as well as the INTEGRAL excess 
within a common framework relying on a single value of S. 

These potential features are interesting, and highlight both the challenges in designing 
a workable cosmology, as well as motivating explorations of some truly novel phenomena. 
While straw man constructions are easy to construct piece by piece, a compelling, unifled 
framework would be worth pursuing. We leave a detailed exploration of these ideas for future 
work. 
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